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X-ray emission from expanding cocoons 
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Abstract. X-ray observations of extragalactic radiosources show strong evidences of interaction between the radio 
emitting plasma and the X-ray emitting ambient gas. In this paper we perform a detailed study of this interaction 
by numerical simulations. We study the propagation of an axisymmetric supersonic jet in an isothermal King 
atmosphere and we analyze the evolution of the resulting X-ray properties and their dependence on the jet 
physical parameters. We show the existence of two distinct and observationally different regimes of interaction, 
with strong and weak shocks. In the first case shells of enhanced X-ray emission are to be expected, while in the 
second case we expect deficit of X-ray emission coincident with the cocoon. By a comparison between analytical 
models and the results of our numerical simulations, we discuss the dependence of the transition between these 
two regimes on the jet parameters and we find that the mean controlling quantity results to be the jet kinetic 
power. We then discuss how the observed jets can be used to constrain the jet properties. 
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1. Introduction 

X-ray observations of extragalactic radio sources have re- 
vealed strong evidences of interaction between the radio 
emitting plasma and the X-ray emitting gas in the ambient 
medium. The observation of Cygnus A with the ROSAT 
HRI (High Resolution Imager) by Carilli, Perley & Harris 
(1994) and with the Chandra X-ray Observatory (Smith 
et al. 2002) showed deficits of X-ray emission in the cluster 
gas spatially coincident with the radio lobes. Observations 
of the Perseus cluster by Bohringer et al. (1993), also with 
ROSAT, and by Fabian et al. (2000) with the Chandra X- 
ray Observatory showed that the cluster emitting gas was 
displaced by the radio lobes of the source NGC 1275. A 
similar behavior has been observed in the Hydra A clus- 
ter, hosting the radio source 3C 218, with the Chandra 
X-ray Observatory by McNamara et al. (2000) (see also 
Nulsen et al. 2002), and in Abell 2052 that shows regions 
devoid of X-ray emission coincident with the radio lobes of 
3C 217 (Blanton et al. 2001). Other clear examples of the 
interaction between radio lobes and the surrounding clus- 
ter gas are given by A 4059 (Heinz et al. 2002, Chandra) 
and A 2199 (Owen & Eilek 1998, ROSAT). Disturbances 
by a radio source are also found in the gas halo of some gi- 
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ant elliptical galaxies such as M 87 (Bohringer et al. 1995, 
ROSAT) and M 84 (Finoguenov & Jones 2001, Chandra). 

On the other hand, theoretical models predict that 
jets in radiogalaxies inflate overpressured cocoons that 
displace and compress the ambient gas and the effects 
of such interaction could indeed expected to be the for- 
mation of cavities and shells in the X-ray emission, as 
shown by observations. A simple one-dimensional model 
of this interaction has been presented by Begelman & 
Ciofh (1989). More detailed and realistic models require 
the use of numerical simulations. In this context, Clarke, 
Harris & Carilli (1997) carried out calculations of the jet 
propagation in a King atmosphere obtaining simulated 
X-ray images to compare with ROSAT data on Cygnus 
A. They demonstrated that a deficit in the X-ray bright- 
ness is indeed shown in the simulation results and found 
agreement between simulations and observations for mod- 
erate Mach number of the jet (M >4). A similar scenario 
is depicted by the numerical simulations of Rizza et al. 
(2000) that showed the interaction and disruption of a 
jet inside a cooling flow cluster. More recently Reynolds, 
Heinz & Begelman (2001) have pointed out that cocoons 
start being strongly overpressured, but, during their evo- 
lution, their pressure decreases, and they then become 
essentially in pressure equilibrium with the ambient or 
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even underpressured. During this evolution, therefore, the 
shock driven in the external medium is strong at the begin- 
ning and becomes very weak at the end. Reynolds, Heinz 
& Begelman (2001) call this last phase "sonic boom". 
The need for weak shocks comes from the observations 
of cool rims surrounding some of the X-ray cavities (see 
A 2052, Hydra A, Perseus A) that rules out the possibil- 
ity of strong shocks driven by the expanding cocoon. To 
explain these observations several analytical (Churazov 
et al. 2000, Soker et al. 2002) and numerical (Churazov 
et al. 2001, Brighenti & Mathews 2002, Quilis, Bower & 
Balogh 2001) models of "bubbles" of hot plasma expand- 
ing subsonically in the ambient medium have been stud- 
ied. If these bubbles are buoyant they can also explain 
the presence of deficits of X-ray emission far from the ra- 
dio lobes as observed in Perseus A. On the other hand, 
Heinz, Reynolds & Begelman (1998) and Alexander (2002) 
have studied self-similar solutions of simplified analytical 
models of overpressured cocoons expanding in a stratified 
medium in order to explain the observed features. 

In this paper we analyze in detail, by using numeri- 
cal simulations, the evolution of the X-ray properties of 
expanding cocoons and their dependence on the jet prop- 
erties. The jets are characterized by their Mach number 
and their density ratio with the ambient medium density; 
the parameter plane is widely covered in order to con- 
sider a wide range of jet powers. We confirm the results 
presented by Reynolds, Heinz & Begelman (2001) on the 
existence of two distinct and observationally different sub- 
sequent regimes of interaction, with strong or weak shocks 
but we are able to determine how and when the transition 
between these two regimes occurs, depending on the jet 
parameters. These results, on the other hand, show how 
the X-ray properties of cocoons could possibly be used as 
diagnostic for the jet characteristics. 

The plan of the paper is the following: in Section 2 
we describe the model, the basic equation and the initial 
conditions, in Section 3 we discuss the X-ray morphologies 
resulting from the simulations, in Section 4 we discuss 
the heating of the external material compressed by the 
expanding cocoon and the consequent changes in the X- 
ray emission properties, in Section 4 we discuss the physics 
of the cocoon expansion that leads to the interpretation 
of the different X-ray morphologies and in Section 5 we 
discuss the astrophysical relevance of our results, finally a 
summary is presented in Section 6. 



stratified (but not self-gravitating) medium 



2. Numerical simulations 

We solve numerically the hydrodynamic equations for a 
supersonic jet, in cylindrical (axial) symmetry in the co- 
ordinates (r, z) continuously injected into a gravitationally 
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where the fluid variables p, p, v and E are, as cus- 
tomary, pressure, density, velocity, and thermal energy 
(p/ (7 — 1)) respectively; 7 is the ratio of the specific heats. 
Radiative losses are neglected since the estimated radia- 
tive times are much longer than the evolution time scale. 

The system of equations (QJ has been solved numer- 
ically employing a PPM (Piecewise Parabolic Method) 
hydrocode (Woodward & Colella 1984). The integration 
domain has a size r domain x z d omain where r domain = 
■^domain = 2.6a, where a is the core radius, defined be- 
low, and has been divided in 1024 x 1024 grid points. The 
axis of the jet is along the left boundary of the domain 
(7' = 0), where we have imposed symmetric boundary 
conditions for p, p, v z and and antisymmetric conditions 
for v r . Reflective boundary conditions are also imposed 
on the boundary of injection of the jet (z = 0) outside its 
radius in order to reproduce a bipolar flow and to avoid 
spurious inflow effects. Free outflow is set on the remaining 
boundaries by imposing a null gradient for each variable 
(d/dr = 0). 

The undisturbed ambient medium is assumed strati- 
fied in a spherically symmetric gravitational well, accord- 
ing to a classical isothermal King model: 
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with R = \/r 2 + z 2 and (3 — 0.5. The resulting pressure 
stratification is kept in equilibrium by an appropriated 
external gravitational potential. 

A (cylindrical) jet is injected from the bottom bound- 
ary of the integration domain, in pressure balance with 
the ambient. The initial jet velocity profile has the form 
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and the corresponding density profile is 
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with w — 0.77, m — 8 and n = 2m, giving a jet radius 
rj = a/20. 

Measuring lengths in units of the core radius a, ve- 
locities in units of the adiabatic sound speed c se in the 
undisturbed external medium and the density in units 
of the ambient central density p(0), our main parameters 
are the Mach number M = v-Jc sc and the density ratio 
v = pj/ p(0). Consistently the unit for the kinetic power is 



£k = ^P(0)r- 
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Table 1. Parameters of the seven simulations performed: 
Mach number M, density ratio v = pj/p(Q), kinetic jet 
power Lj/Lk = M^v and final time of the simulation r^n = 

tfin/to 



M 


V 


Lj/Lk 


Tfln 


10 


0.1 


10 A 


1.38 


60 


0.001 


2.16 x 10^ 


1.22 


60 


0.01 


2.16 x 10^ 


0.49 


60 


0.1 


2.16 x 10 4 


0.15 


120 


0.001 


1.73 x 10^ 


0.76 


120 


0.01 


1.73 x 10 4 


0.25 


120 


0.1 


1.73 x 10 b 


0.07 



= L5xl ° 4 1ilF^) (50l^) 2 (^v) 3/2ergs " 

(5) 

and the unit of time is 

*° = ^ 8 * 10? fe) < 6 > 

where uq is the central electron density, and we will mea- 
sure time t in units of tg. Then the kinetic power of the 
jet expressed in unity of is given by 

Lj = M 3 j/L k . (7) 

In Tabled we show the parameters M and v for the simu- 
lations that we have performed giving also the correspond- 
ing values of the kinetic jet power Lj in units of Lk and 
the final simulation times in units of to. 



3. X-ray signatures: shell and cavity 

The general structure of the interaction between a low 
density jet and the ambient medium is well known since 
the first simulations of Norman et al. (1982) (see also 
Massaglia, Bodo & Ferrari 1996 and Krause 2003) ). The 
flowing jet matter, slowed down by one or more termi- 
nal shocks, inflates a cocoon that compresses and drives 
shocks in the surrounding external medium. The com- 
pressed ambient material forms a shell surrounding the 
cocoon: the boundary between the shell and the cocoon 
is marked by a contact discontinuity, while the boundary 
between the shell and the undisturbed external medium 
is marked by a shock. We will call 11 extended cocoon" the 
whole region interested in the interaction between the jet 
and the ambient medium. The ll extended cocoon" is then 
formed by the cocoon proper and by the surrounding shell. 
The cocoon, which is formed by the expanded jet material 
forms a cavity with very low density and high tempera- 
ture, in which the X-ray emissivity is strongly depressed. 
On the other hand the external material in the shell has an 



enhanced emission due to its compression. In addition, de- 
pending on the external shock strength this material can 
be heated and its emission properties can then change. 

In principle, we then expect three main features in the 
X-ray properties of the region of interaction between a jet 
and the ambient medium: 

1. a region of depressed emission coincident with the co- 
coon; 

2. a shell of enhanced emission; 

3. a variation in the spectral properties of the emission 
by the shell. 

The actual appearance, however, will depend on an inter- 
play between all these effects. In fact, we have to take into 
account that what we see is an integration along the line 
of sight, that can go both through the cocoon and through 
the shell. Therefore, it can happen that the enhanced shell 
emission can compensate the deficit of emission in the co- 
coon and, in these cases, in correspondence of the cocoon 
we can even see a flux higher than that expected in the 
undisturbed case. 

In Figs. (^-13) we show in the first two columns simu- 
lated X-ray flux distribution for the seven cases described 
in Table ^ The figures are symmetrical with respect to 
the r and z axes. The emissivity per frequency unit is 
computed with a Raymond-Smith thermal spectrum code 
(see HEASARC Web Page) from the electron density and 
temperature distribution obtained in our calculation. The 
emissivity e = n 2 A(T) is then integrated in the 0.1 — 4 keV 
band. The flux distribution is calculated integrating the 
optically thin emissivity along the line of sight A that is 
assumed perpendicular to the jet axis: 

/+rdomain 
e(n,T,\)d\. 
-I-domain 

We consider an isothermal ambient medium with a tem- 
perature T = 2.3 keV and a central electron density 
no = 0.01 cm -3 . In FigQ]we show the results for the three 
cases with M = 120, each row refer to a different value 
of the density ratio and more precisely the top row is for 
v = 0.1, the middle row is for v — 0.01 and the bottom 
row is for v — 0.001. The first two panels in a row are for 
a different time and we have chosen the times in order to 
have lengths respectively of 1 and 2 core radii. Figure El 
has the same structure, but is for the cases with M = 60, 
while Fig-EJis f° r the single case with M = 10 and v = 0.1. 
A more quantitative view of these results can be obtained 
presenting cuts of these figures along selected directions. 
In particular in Figs. Q-|3J in the third and fourth col- 
umn we show longitudinal cuts of the images along the jet 
axis. The figures correspond to the images of the first two 
columns and refer respectively to M — 120, M — 60 and 
M — 10. The dashed curves in all the figures represent 
the emission from the undisturbed atmosphere. 

These figures show a sequence of morphologies starting 
from cases in which a cavity is not present going to cases 
in which the cavity is the dominant feature of the image. 
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Fig. 1. Simulated X-ray fluxes in the 0.1 — 4 keV band for the M = 120 cases. The rows refer from top to bottom 
to the v = 0.1,0.01,0.001 cases respectively. In the first two columns on the left the X-ray images are shown in 
logarithmic scale at times corresponding to cocoon lengths of l c ~ 1 and 2 core radii. In the two columns on the right 
longitudinal cuts along to the jet axis corresponding to the images on the left are plotted with a solid line. The dashed 
line represents the emission from the undisturbed atmosphere. The fluxes are given in unity of the central flux of the 
unperturbed atmosphere 



In particular, we see that for M — 120 the brightness 
depression is evident only for the lighter case, while for the 
M — 60 cases starts to be present for v — 0.01 during its 
evolution and is very evident during the whole evolution 
for v = 0.001. For M = 10 the cavity is dominant already 
for v = 0.1. 



The presence or absence of the brightness depression, 
as discussed by Clarke, Harris & Carilli (1997) depends on 
the thickness of the shell. In fact, the line of sight that goes 
through the cocoon region with very low emissivity crosses 
also the shell of enhanced emissivity, and the observed flux 
is higher or lower than the undisturbed profile depending 
on the interplay between the two effects. If the shell is 
narrow, as a consequence of mass conservation, it will have 
a high density, the emissivity will be greatly enhanced 
and will overcome the decrease in the cocoon giving an 
observed brightness higher than the undisturbed profile. 
Following Clarke, Harris & Carilli (1997), we can write 



the ratio of the observed flux to the undisturbed one as 



/ 5(2- 5f ' 

where we have neglected the dependence of emissivity on 
the temperature and 8 is the ratio between the shell width 
and the cocoon radius. From Eq. JSJ), we see that the ratio 
increases as 5 goes to and becomes less than 1 for 5 > 
0.38. 

In Fig. 0] we show the density distribution for the 
M = 60 cases: each row refer to a different density ra- 
tio (y = 0.1, 0.01, 0.001 from top to bottom) and each col- 
umn is for different times corresponding to cocoon lengths 
1,1.5,2 respectively. We see in fact that the width of the 
shell increases going from the high to the low v case. 
Moreover we can see that the shell tends to widen and 
the bow shock becomes weaker during the evolution of 
the cocoon, a part from the v = 0.1 case in which the rel- 
ative thickness of the shell remains constant. Note that the 
localized feature appearing along the z — axis is due to 
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Fig. 2. Simulated X-ray fluxes, in logarithmic scale, in the 0.1 — 4 keV band for the M = 60 cases. The rows refer 
from top to bottom to the v = 0.1, 0.01, 0.001 cases respectively. The quantities represented are the same as in Fig.^ 
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Fig. 3. Simulated X-ray fluxes, in logarithmic scale, in the 0.1 — 4 keV band for the M = 10, v = 0.1 case. The 
quantities represented are the same as in Fig. 1 1121 



the assumed reflection boundary conditions and does not 
affect the general morphological and dynamical behavior. 

In order to better quantify the deficit or the enhance- 
ment in the X-ray emission, we have computed integral 
measures K c and K s of these quantities defined as 

Ko = Jr J (/ - h)dS (9) 
K B = 4" / (/ - fk)dS , (10) 



where the surface integrals are computed over the domains 
S c and S'g defined as the areas over which the integrand 
(/ — /k) is respectively lower (the cavity) and greater (the 
shell) than zero, where / is the flux at a given position and 
/k is the flux at the same position for an unperturbed King 
atmosphere. The quantities K c and K s are calculated in 
two different energy bands (0.1—4 keV and 4—10 keV) and 
are plotted as a function of the length of the cocoon l c (t) in 
Fig. El the upper panels refer to K c , while the lower panels 
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Fig. 4. Density maps for the M — 60 cases. The scale is logarithmic. The rows refer from top to bottom to the 
v = 0.1,0.01,0.001 cases respectively. In the columns density maps are shown at times corresponding to a cocoon 
length l c ~ 1, 1.5 and 2 core radii. The densities are given in unity of the central density of the unperturbed atmosphere 



refer to K s . In the figure we present the results for all 
the values of v considered and for a single representative 
value of the Mach number M = 60. In each figure the 
solid curve is for v = 0.1, the dashed curve is for v = 0.01 
and the dashed-dotted curve for v — 0.001. Looking at 
the figure for K c (0.1 — 4keV), we see that, as expected 
from the results discussed above, in the case of v = 0.1 
the deficit is very low and does not increase with time, 
instead, the deficit increases at the decreasing of v and 



shows an increase with time. When we look at the excess 
emission in the 0.1 — 4 keV band, measured by K Sl we see 
that the cases v — 0.001 and v = 0.01 behave as expected, 
with the case v — 0.01 presenting an higher shell emission 
with values decreasing in time (remember that the shell 
widens and decreases its density during its evolution) . The 
case v = 0.1, instead, presents an unexpected behavior. 
In fact its value of K s is lower than that found for v = 
0.01 and stays almost constant with time. This however 
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Fig. 5. (Upper panels) Plot of the integrated flux deficit 
as a function of cocoon length in the 0.1 — 4 keV band 
(left) and in the 4-10 keV band (right) for the M = 60 
cases. This quantity is defined in Eq. 10. (Lower panels) 
Plot of the excess emission as a function of cocoon length 
in the 0.1 - 4 keV band (left) and in the 4-10 keV band 
(right) for the M — 60 cases. This quantity is defined in 
Eq. I |10| l. In the four panels the three lines refer to the 
v = 0.1 (solid), v = 0.01 (dashed) and v = 0.001 (dash- 
dotted) cases 

can be understood by considering that the shock in the 
external medium heats this medium changing its emission 
properties. The spectral range we are considering for the 
flux calculation is 0.1 — 4 keV and the material in the 
shell may have been heated to a temperature for which the 
maximum of the emission falls in a harder spectral range. 
These effects will be however analyzed in more detail in 
the next Section 0] 

4. The shell temperature 

As we discussed, the shock driven in the external medium 
can heat it and its effects can be more or less evident 
depending on its strength. The heating of the shell can 
change the typical emission energies, an example of the 
consequences of this effect has been shown in the bottom 
panel of the first column of Fig. [3] where we have seen that 
the case with v = 0.1 seemed to present an anomalous be- 
havior. We can now compare the first column of Fig. 
with the second column where we show the same quanti- 
ties but in the range 4—10 keV. Looking at the bottom 
panels we see that the case v — 0.1, which in the softer 
band presented an excess emission below that of the case 
v = 0.01, shows in the harder band an excess emission 
larger than the other cases. The gas in the compressed 
shell becomes in fact hotter as we increase the value of 
v and the emission is then shifted towards higher emis- 
sion energies. To quantify in more detail this temperature 
change, we have plotted in the first row of Figures EE the 
quantity 



where ?i(x) is the temperature at a given position, 5(T\ — 
T) is the Dirac delta function and e are the radiative losses 
per unit volume in the 0.1 — 10 keV band as calculated 
with the Raymond-Smith code. This quantity measures 
the energy emitted by the gas at the temperature T in 
the 0.1 — 10 keV band. In in the first row of Fig. 00 and 
[S]we have plotted this quantity for the Mach numbers 120, 
60 and 10 respectively. In each figure columns are for dif- 
ferent values of v and the panels in each column are for 
the time corresponding to a cocoon length of 2. In each 



E(T)= f e(n,T 1 ,x)5(T 1 (x) 
Jv 



0.07 



0.69 



1 10 
T (Kev) 



1 10 
T (Kev) 



1 10 

T (Kev) 





20 40 60 80 
dee 



20 40 60 80 
deg 




(11) 



Fig. 6. Temperature plots for the M — 120 cases (v = 
0.1,0.01,0.001 from left to right respectively) at a time 
corresponding to a cocoon length of 2 core radii. ( Upper 
panels) Plot of the quantity E(T) as defined in Eq. (fT/T/f) . 
It represents the energy emitted in the 0.1 — 10 keV band 
by the gas at the temperature T. The vertical line marks 
the temperature of the ambient isothermal medium while 
the peak correponds to the compressed shell emission. The 
emission is normalized to the total emission from the un- 
perturbed atmosphere. (Lower panels) Plot of the shell 
average temperature as a function of the angle with origin 
in r = 0, z = 0. The 0° angle corresponds to the direction 
of the jet axis. The ambient temperature is taken 2.3 keV 



panel we see a vertical line marking the temperature of the 
ambient isothermal medium and a peak correponding to 
the emission from the shocked shell material. We can see 
that the temperature of the shell increases as we increase 
the Mach number and the density ratio v. Taking into ac- 
count the cases with the lowest density ratio v = 0.001 we 
can see that in the M — 120 case the shell temperature is 
a factor two higher than the ambient one, in the M = 60 
case it is slightly higher while in the M = 10 case it is 
clear that there is emission also from gas with a temper- 
ature lower than the ambient one. This fact agrees with 
several observations (e.g. Perseus A) in which the coolest 
gas is found in the shell of enhanced emission. The cool- 
ing in this simulation is due to the adiabatic expansion 
of the shell after being compressed by a weak shock. As 
time elapses, radiative losses may become important and 
contribute to the cooling of the shell: for a central density 
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Fig. 7. Temperature plots for the M = 60 cases [y = 
0.1,0.01,0.001 from left to right respectively) at a time 
corresponding to a cocoon length of 2 core radii. The quan- 
tities plotted are the same as in Fig. El 
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Fig. 8. Temperature plots for the M = 10, v = 0.1 case 
at a time corresponding to a cocoon length of 2 core radii. 
The quantities plotted are the same as in Fig. HJ1 and in 
Pig. El 



of 0.04 cm' 3 and a temperature of 3keV (e.g. Perseus A, 
see Schmidt, Fabian & Sanders 2002), the radiative cool- 
ing time is ~ 3 x 10 7 years, i.e. close to the simulated time 
of the M = 10, v — 0.1 case. The shell temperature can 
reach temperature up to 800 keV for the case M = 120 
and v = 0.1. The heating of the shell is not uniform but 
depends on the shock strength and it is higher towards the 
jet head and becomes lower further from the head. This is 
represented in the second row of Figs. 10181 where we show 
the average shell temperature as a function of angle for all 
the cases considered. 

From the figures we see that, for the cases with v = 0.1 
and M — 60, 120, the shell temperature is everywhere 
larger than 8 keV and reaches temperatures up to 850 keV 
at the jet head of the M = 120 case. The lower density 
cases, instead, present an increase in temperature mainly 
concentrated in the forward part of the cocoon. This tem- 
perature distribution will have consequences in the shell 



morphology as seen in different X-ray bands. In fact, we 
expect to see emission from the forward part of the cocoon 
at high energies, while the backward part is expected to 
be more dominant at lower energies. 

5. The cocoon dynamics 

The numerical simulations presented in previous sections 
have shown that the expanding cocoon drives a shock 
in the ambient medium, which is then compressed and 
heated. The amount of compression and heating is fun- 
damental in determining the X-ray emission properties. 
The results we have shown tell us that, as also discussed 
by Reynolds, Heinz & Begelman (2001), the strength of 
the shock, driven by the cocoon, weakens as the cocoon 
expands and the stage at which the transition between a 
strong and a weak shock occurs depends on the jet prop- 
erties, i.e. its Mach number and its density ratio. In this 
section we will try to determine in a more quantitative 
way this dependence of the transition on the jet physical 
properties. We will do that through a comparison between 
analytical models for the cocoon expansion in an homo- 
geneous medium and the results of our numerical simula- 
tions. In the next subsection we will then examine the co- 
coon dynamics in a uniform medium, in subsection 15 . 21 we 
will examine the dynamics in a stratified medium and in 
in subsection 15 . 31 we will try to determine the dependence 
of the transition discussed above on the jet properties. 

5.1. Expansion in a uniform medium 

The first attempt to build an analytical model for the co- 
coon dynamics is due to Begelman & Cioffi (1989) and 
Cioffi & Blondin (1992). They considered only the case of 
strongly overpressured cocoons, for which they consider 
only strong shocks driven in the external medium and they 
essentially neglect the external pressure. For our purposes, 
we have to extend the description to a more general situ- 
ation, in which the external shock can be of any strength 
and the external pressure is taken into account. We de- 
scribe in detail the model in the Appendix, here we give 
only the resulting behavior of the extended cocoon ra- 
dius and cocoon pressure versus time: the cocoon radius 
is given by 



r 2 = c 2 se t 2 + vV - lrX /4 Mc se i 



(12) 



and the cocoon pressure is 



7(7- 1) 



1 



„2 f 2 



1 C se £ 



U 1/2 M 2 r 2 



v 1 / 4 M n 



(13) 

Equations 1121 1131) describe two different phases in the 
evolution: initially the cocoon is strongly overpressured 
and the solutions behave in the same way as described 
by the Begelman & Cioffi (1989) model, with the radius 
proportional to t 1 ^ 2 and the pressure proportional to t -1 . 
As the pressure decreases, the contribution of the external 
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pressure becomes more important and the external shock 
becomes weak. In this second phase the pressure tends to 
become constant and the extended cocoon expands essen- 
tially at the sound speed and increases as cx t. We can use 
Eq. I|A.4() to determine the radius of the cocoon proper: in 
the first phase the internal radius increase cx t l/2 like the 
external one not allowing the shell to expand, while in the 
second phase it tends to a constant value. The compression 
of the external medium can be related to the relative shell 
thickness and we can observe that, while it stays constant 
in the first phase, as already discussed above, it tends to 
increase in the second phase. The transition between these 
two phases happens when the two terms on the right hand 
side of Eq. (|T2|) become comparable; then, assuming a co- 
coon length / c = v l / 2 Vtt, we can determine a scaling law 
for the length of transition I* between the two regimes 



I* cx MV/ 4 



(14) 



Therefore, for a given Mach number and for a given length 
of the cocoon, we expect the jets with lower density ratios 
to be less overpressured and therefore to form wider and 
less dense shells. 

In a similar way we can also find a solution assum- 
ing a spherical symmetry for the cocoon, that is solving 
the system of equations (|A.6I - lA~%)) taking V c — 47r/3r 3 . 
Taking into account separately the strongly overpressured 
(P c ^> P*) and the weakly overpressured (P c ~ P*) 
phases we can find two distinct regimes of expansion: a 
supersonic one during which the external radius behaves 
like 



25( 7 2 - 1) 
9?r 



1/5 



{- 

\Poxt 



1/5 



f 3/5 



and the pressure decreases as 



Pc 



3(7-1) 

47T 



9 71 



25( 7 2 - 1) 



3/5 



(A 



4/5 



(15) 



(16) 



and a weakly overpressured phase during which the cocoon 
pressure tends to a constant value and the radius r c ex- 
pands cx t at the sound speed. Assuming spherical symme- 
try also for the cocoon proper (V C i — 47i3r?) we can derive 
the behavior of the cocoon radius r; solving Eq. (|A.5(I : rj 
expands as i 3 / 5 during the strongly overpressured phase, 
with the same rate of expansion of r c , while it behaves as 
i 1 / 3 during the weakly overpressured one. As in the cylin- 
drical geometry model the shell is allowed to widen only in 
the weakly overpressured phase. We can then determine 
a scaling law for the radius of transition between the two 
regimes imposing that the expansion speed determined by 
Eq. (|15JI approaches the sound speed c se 



C soPoxt 



1/2 



(17) 



In the spherical symmetry model the cocoon properties 
scale with the jet power since it is the only jet parameter 
that enters the system of equation solved in the Appendix. 



It is worthwhile noticing that our solution for the spher- 
ical strongly overpressured cocoon (Eq. ^] and 1161) is the 
same, apart from numerical constants, of Heinz, Reynolds 
& Begelman (1998) who solved a system of equations sim- 
ilar to ours but limited their solution to the strongly over- 
pressured regime. 

Looking at the geometry and at the detailed structure 
of the shell and the cocoon, these analytical models are 
too simplified to describe in detail the radiative proper- 
ties of the cocoon but still they show clearly how the rel- 
ative thickness of the shell can increase only in a weakly 
overpressured regime. 

5.2. Expansion in a stratified medium 

The cocoon expansion in a stratified medium is, of course, 
much more complicated, since the external pressure is not 
constant, and the development of an analytical model be- 
comes more problematic. However the above discussion 
can provide a framework for understanding the results of 
numerical simulation. 

We recall that the presence of a brightness depression 
depends substantially on the relative thickness of the shell. 
We then must first study the evolution of the extended co- 
coon radius and the cocoon proper one. We first analyze 
the behavior of the extended cocoon radius r e by compar- 
ing the results directly obtained by the numerical simula- 
tions with an estimate obtained, using the average cocoon 
pressure as the driver for the expansion. More precisely, we 
obtain this estimate by integrating numerically Eq. (|A.8|) , 
where for P c (t) we use the average value obtained at every 
time step from the simulations. For defining the extended 
cocoon radius in the simulations, we concentrate on the 
base portion of the cocoon and we take an average value 
of the bow shock radius between 1/8 and 1/4 of the total 
length of the cocoon. The results of the comparison are 
represented in Fig. 03 where we plot r and r; as a func- 
tion of time for all the cases we have considered. In all the 
panels the solid curves correspond to the results obtained 
in the simulations, while the dashed curves correspond to 
the estimates obtained through the average pressure. The 
figures show that this estimate reproduces very well the 
actual behavior of the extended cocoon radius: this result 
tells us that the average cocoon pressure is a good esti- 
mate of the local pressure driving the cocoon expansion. 
In a similar way we can proceed for the cocoon radius n, 
for which we can use eq. I|A.5(I where with rjLj we intend 
the fraction of jet power that is thermalized and goes into 
internal energy of the cocoon. This fraction can be esti- 
mated from eq. (|A.6|I . Eq. (|A.4|I can be integrated giving 



7 c = - 7 i I77 ^P c (0) 1 /7 Fc (o) + l_i / tjL^-^dt 



Po(t) 



7 Jt 



(18) 

where again for P c we can use the average cocoon pres- 
sure obtained above. This equation gives us an estimated 
behavior of the cocoon volume and from this we can esti- 
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mate the behavior of the cocoon radius as r\ = \JV c jKl c 
where K is a form factor that depends on the shape of 
the cocoon. One problem in determining r- x is that the 
shape of the cocoon is strongly variable since the interface 
between cocoon and external shocked material is subject 
to Kelvin-Helmholtz instabilities. Moreover, the shape of 
the cocoon may depend on the jet physical parameters. 
However, we can see from Fig. (where, as before with r e , 
we compare this estimate with the actual results obtained 
from the simulations) that using an appropriated average 
values for K, we can capture quite well the average behav- 
ior. The form factor should be K = 7r for a cylinder and 
K = 2ir/3 for a hemisphere. The value that can be esti- 
mated from the simulations and that works quite well for 
every case is K ~ 2. Again, this good agreement between 
the estimates obtained through the usage of the average 
cocoon pressure and the actual results tells us that we can 
use the average pressure behavior for understanding the 
cocoon dynamics 

Looking in more detail at the behavior of the extended 
cocoon radius, we observe that, as expected from the dis- 
cussion of the uniform case, its expansion velocity de- 
creases faster initially, when the cocoon is strongly over- 
pressured, and then decreases more slowly, when the shock 
becomes weaker. Fitting a power law oc r Q for the quan- 
tities l c , r c , r*i in the initial and terminal part of the evo- 
lution for the cases with M = 60, we get the exponents 
a reported in Table We see that, at later times, the 
exponents of r c become systematically larger than those 
obtained at the initial times. Moreover we see that the 
differences between the two exponents are larger for the 
lower density cases and for v — 0.1 are minimal. In this 
last case, in fact the cocoon stays strongly overpressured 
during the whole evolution. 



Table 2. The Table shows the exponent a of the fitted 
power laws r Q for the quantities l c (cocoon length), r c 
(extended cocoon radius) and r; (internal cocoon radius) . 
The left column refers to the initial phase of the evolu- 
tion of the cocoon (l c < 1.0) while the second one to the 
advanced one (l c > 1.5) 
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Fig. 9. Plot of the internal and external radii vs time. 
The rows refer to the different Mach cases while columns 
to the different v values. The solid lines show the extended 
and internal cocoon radii as determined in the simulations 
taking an average value between 1 /8 and 1 /4 of the total 
length of the cocoon. The dashed lines show the same 
quantities as determined by our analytical model assuming 
that the expansion of the cocoon is driven by the average 
pressure of the cocoon 

A similar analysis can be done for the cocoon radius r\. 
We see that in the high and intermediate density cases the 
exponent for r\ also increases but less than the exponent 
for r e . In the low density case, instead, T\ becomes almost 
constant in time. This behavior tells us that, when the 
cocoon is not any more strongly overpressured, the rela- 
tive shell thickness start to increase as it happened in the 
homogeneous case. The behavior of the relative thickness 
of the shell can be in fact derived from Table J2J), since 
it is related to the ratio r-Jr c . We then see that it has in 
general a slower increase at the beginning of the evolu- 
tion, and accelerates in the following phases. As discussed 
above, the effects of not being strongly overpressured are 
more evident by decreasing the value of v and, in fact, the 
increase of the relative cocoon width is larger for smaller 
v. If we compare the cocoon widths for equal distances of 
jet propagation we have a further effect that amplify the 
consequences of the behavior discussed above, namely we 
have to take into account that the advance velocity of the 
head of the jet decreases when we decrease the value of v. 

In this subsection we have seen that the transition 
from a strongly overpressured cocoon to a cocoon which is 
essentially in pressure equilibrium with the surroundings 
leads to different properties of the shell of the compressed 
external material. In the next subsection we will try to 
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see whether it is possible to determine in a more quanti- 
tative way how this transition depends on the jet physical 
parameters. 

5.3. The transition from strong to weak shocks 

From the above discussion we have seen that the aver- 
age cocoon pressure can be used for interpreting the co- 
coon dynamics and therefore determining the transition 
between the strong and weak shock regimes. In the case 
of uniform external medium, the cocoon pressure, after a 
decrease proportional to t , tends to a constant. In the 
present case, with a decreasing external pressure, we ex- 
pect that the cocoon pressure does not tend to a constant 
but that it will tend to follow the behavior of P* , defined 
in eq. I|A.7|) . 
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Fig. 10. Plot of the pressure vs cocoon length. The rows 
refer to the different Mach numbers while columns to the 
different v values. The solid line shows the evolution of the 
average cocoon pressure with time in the different cases. 
The dashed one represents the external pressure averaged 
on the cocoon volume as defined in Eq. I A. 71 



strongly overpressured to weakly overpressured cocoons 
for M = 120, v = 0.001 and M = 60, v = 0.01; fi- 
nally we have two cases where the cocoon is always weakly 
overpressured, i.e. for M — 60, v = 0.001 and M = 10, 
v = 0.1. This is in agreement with the discussion on the 
X-ray flux distributions following Figs. Q-|3J). 

In order to determine the length of the cocoon at which 
the transition occurs, we look for a scaling law for the 
quantity (P c — P*)/P* as a function of the parameters l c , 
M and v during the initial strongly overpressured phase. 
Taking into account only the cases that are stronlgy over- 
pressured at the beginning of the evolution, we notice that 
the initial decrease of the average cocoon pressure follows 
a similar behavior for all the cases that we have consid- 
ered: with a power law fit (P c - P*)/P* oc l~ a to the ini- 
tial evolution of the different cases we find a mean value 
a = 0.9 with an uncertainness of 10%. In Table 01 we 
show the values of the quantity (P c — P*)/P* at the be- 
ginning of the evolution, when the cocoon length is one 
third the core radius and the cocoon is typically strongly 
overpressured (notice that the cases M = 10, v = 0.1 and 
M = 60, v = 0.001 make an exception). These values can 
be represented, with good approximation by the scaling 
(P c - P*)/P* oc M 1 ' 85 ^ - 62 . The general scaling law thus 
becomes: 

Pc ~ P * oc /-o.9 M i.88 I/ o.ra _ (19) 

The scaling relation (|19fl shows that the initial pressure of 
the cocoon is not proportional to an arbitrary combina- 
tion of the jet parameters M and v but it scales with the 
normalized kinetic jet power L-JL^ = M^v (see Eq. Q) 
as (P c - P*)/P* oc l-°- 9 (Lj/L k ) - 62 . This result is con- 
firmed by the plot of the values of (P c — P*)/P* reported 
in Table against the jet kinetic luminosity Lj/L^, after 
setting l c — 0.3. The plot is shown in Fig. 1111 

Table 3. Values of the quantity (P c - P*)/P* for the 
different simulations evaluated at a cocoon length l c = 0.3 
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In Fig. E| we have plotted the behavior of P c (solid 
curve) and P* (dashed curve) against the cocoon length, 
where P c is the average cocoon pressure and P* is the 
external pressure averaged over the cocoon volume. The 
ratio between the two quantities gives a measure of the av- 
erage strength of the shock driven in the external medium. 
Looking at the behavior of the cocoon evolution in the 
different cases we have an always strongly overpressured 
expansion for the cases M — 120 and v — 0.1, 0.01 and 
for M = 60 and v — 0.1, in these three cases the transi- 
tion length would be larger than the actual longitudinal 
size of our domain; two cases showing a transition from 



Of course the transition from a strongly to a weakly 
overpressured regime is not sudden, but gradually happens 
when (P c — P*)/P* becomes of the order of unity. Setting 
a constant value for (P c — P*)/P* of this order in Eq. lj*H?|l 
we derive the scaling law for transition length I* (as we 
have done for the uniform case, see Eq. dJ). 

II oc M 2 - 05 v - 69 (20) 
or in terms of the jet kinetic power 

/ T \ 69 

'Ht) (21) 
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Fig. 11. Plot of the values of (P c - P*)/P* at a cocoon 
length l c = 0.3 (see Table against jet kinetic luminos- 
ity Lj/L k = M 3 v. The solid line is a (P c - P*)/P* cx 



(Lj/Lk) , that is also the best fit, excluding the value 
of the case M = 10, ^ = 0.1 that is in a weakly overpres- 
sured regime already 



Equation (|20J) represents a family of curves in the plane 
(M, v) , and, for a given cocoon length one of the curves 
of this family will mark the separation between cocoons 
that drive strong shocks in the ambient medium (strongly 
overpressured regime) and cocoons that drive only weak 
shocks (weakly overpressured regime). In the next section 
we will discuss the astrophysical relevance of these results. 

6. Discussion 

In the previous sections we have seen that strongly and 
weakly overpressured cocoons present different X-ray mor- 
phologies. The former ones will show no deficit of emission 
accompanied by a strong emission from a shell marking 
the shock, driven by the cocoon expansion. The material 
in the shell will be much hotter than the ambient medium, 
the shell emission will then be shifted to higher frequen- 
cies and therefore more visible at higher X-ray energies, 
especially near the jet head. Weakly overpressured cocoons 
will be instead characterized by the presence of a deficit 
of emission in the cocoon, while the emission from the 
shell will be much less visible than in the previous case; 
in addition, the material in the shell will be essentially at 
the same temperature as the ambient medium and there 
will be no change in the emission spectrum. These two 
regimes depend on the jet physical parameters and on the 
age of the cocoon. High Mach number jets with high den- 
sities will be strongly overpressured for a longer fraction of 
their life, decreasing the density and the Mach number the 
transition to weakly overpressured cocoon will occur at an 
earlier stage. We have found that the scaling of the transi- 
tion length from one regime to the other is given to a good 
approximation by Eq. I|20(l . Fixing a length of the cocoon, 
from this equation we can then obtain a relation between 
Mach number and density ratio such as M = C(? c )^ _1 / 3 . 



Lj/L^lO 4 Lj/L^lO 6 




Fig. 12. Representation of the different regimes of the co- 
coon expansion in the M — v plane. Inside the dashed area 
the jet becomes subsonic with respect to its internal sound 
speed. The solid lines correspond to a constant kinetic jet 
power Lj/Lk = M 3 v given in unit of the power Lk defined 
in Eq. The fragmented curves separate strongly and 
weakly overpressured regimes for different lengths of the 
cocoon given in core radii units. It can be seen that these 
lines correspond to lines of constant kinetic jet power. The 
diamonds correspond to the simulated cases 



This scaling law has been represented in Fig. ^] with the 
dotted, dashed and dash-dotted lines for cocoon lengths 
equal to 0.5, 1 and 2 core radii respectively. The propor- 
tionality constant C(l c ) has been set to fulfill the transi- 
tion conditions in the cases where the transition between 
the two regimes is observed, i.e. M = 60, v — 0.01 and 
M = 120, v = 0.001. In the same figure we can also notice 
the presence of another line, that individuates a region in 
which the jet becomes subsonic with respect to its internal 
sound speed. 

Since Mach number and density ratio are physically 
very meaningful parameters, but cannot be determined 
for actual jets, it is favorable to express the scaling re- 
lation in terms of the jet kinetic power. In order to do 
so, we make use of Eq. (|21|) where Lk is a quantity that 
depends on the properties of the environment. From this 
relation we see that jets with the same kinetic power have 
the same transition length between the two phases. The 
curves in Fig. 1121 separating the strongly and weakly over- 
pressured regimes, which are given by M cx i/ -1 / 3 , will 
therefore correspond also to a constant value of Lj/Lk- 
These considerations tell us that the separation between 
the two regimes is essentially determined by the jet ki- 
netic power: the transition will occur at higher values of 
l c for high kinetic power jets and at lower l c for low kinetic 
power jets. For instance the dashed curve shown in Fig. 
1121 representing the separation between the two regimes 
for a cocoon length equal to 2 core radii corresponds to 
Lj/Lk ~ 4.7 x 10 3 . The properties of the environment en- 
ter in the parameter Lk that fixes a measure for the jet 
kinetic power. For a typical cluster environment with a 
central density 10~ 2 cm~ 3 , a core radius 100 kpc and a 
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temperature 3 keV we have = 6 x 10 42 ergs -1 and the 
dividing kinetic power for a cocoon length equal to 1.5 
core radii will be Lj ~ 1.9 x 10 46 ergs -1 . It is important 
to notice that for jets that are slightly supersonic with re- 
spect to their internal sound speed the jet enthalpy flux 
becomes comparable to the kinetic one. In these cases our 
estimates for the jet kinetic power should be corrected at 
most by a factor two in order to obtain an approximation 
of the total power of the jet. 

Considering the particular case of Cygnus A, its clus- 
ter environment is characterized by a central density 
0.07 cm" 3 , a core radius 35 kpc and a temperature 3.4 
keV (Cariili et al. 1994), giving L k = 6.2 x lO^ergs" 1 . 
The radio lobes of Cygnus A show an extent of ~ 70 kpc, 
twice the core radius (Z c = 2). For this length the dividing 
power is given by L-JL^ ~ 4.7 x 10 3 , as it is possible to see 
in Fig. E| Since Cygnus A clearly shows deficit of X-ray 
emission and therefore is in a weakly overpressured regime, 
these estimates give an upper limit to the kinetic power 
of its jet ~ 2.9 x 10 46 ergs -1 . Nevertheless the pressure 
of the expanding cocoon must be higher than the ambi- 
ent one since observations by Chandra (Smith et al. 2002) 
show clearly that the shell is slightly hotter than the sur- 
rounding medium. Then the cocoon is still expanding as 
a (weak) shock wave more than a sound wave. 

Another well known example of the interaction be- 
tween a radio source and the thermal gas of a cluster is the 
FR I type radio galaxy 3C 84 inside the Perseus cluster. 
The cluster core can be modeled with a central electron 
density 0.04 cm -3 , a core radius a ~ 50 kpc and a central 
temperature T ~ 3.1 keV (Schmidt, Fabian & Sanders 
2002) yielding — 6 x 10 42 ergs -1 . The radio lobes show 
an average extent of ~ 22 kpc, giving a dividing power 
Lj/Lk ~ 5.2 x 10 2 . These estimates give an upper limit on 
the kinetic luminosity of the jet Lj = 3.1 x 10 45 ergs -1 . 
This upper limit can be lowered considering that the shell 
of enhanced emission surrounding the cavities contains 
clearly the X-ray coolest gas in the cluster (Fabian et al. 
2001). This fact rules out the presence of strong shocks 
driven by the expanding cocoon. In our scheme this means 
that the cocoon must be in a weakly overpressured regime 
in the first phase of expansion yet as for example in the 
case M — 10, v = 0.1. As it can be seen in Fig. [S] the av- 
erage temperature of the shell for this simulation is equal 
to the ambient one and in the shell there is also emission 
from gas cooler than the ambient one. This cooling is due 
to the adiabatic expansion of the shell after being com- 
pressed by a weak shock, in a way similar to that described 
in the simulations of expanding hot bubbles by Brighenti 
& Mathews (2002). Taking this case [L-jLy = 100) to de- 
termine a limit on the kinetic power of the jet we obtain 
Lj = 6. 3x 10 44 ergs -1 . For deriving the total jet power this 
estimate should be corrected by the enthalpy term, that 
for the above values of M and v is approximately equal 
to the 30% of the kinetic power. This estimate is similar 
to the one by Fabian et al. (2001) who used the analyti- 
cal bubble model by Churazov et al. (2000) to determine 
their limits. This substantial agreement can be understood 



considering that the Churazov et al. (2000) model solves 
exactly Eq. I|A.4|) with spherical symmetry assuming that 
the bubble pressure is equal to the external one. The limits 
of Fabian et al. (2001) are then determined requiring that 
the bubble has the observed dimensions, that it expands 
subsonically and that it is not buoyant. This situation is 
similar to what we refer to as a weakly or not overpres- 
sured regime. The great difference between our low power 
cases and a hot underdense bubble is the presence of the 
highly collimated jet that forms a thermal hot spot and 
is likely not effective in uplifting the cooling flow gas as 
proposed for example by Soker et al. (2002) to explain 
the presence of cool gas in the rims around the cavities. 
A collimated jet tends to go through the gas of the clus- 
ter displacing it aside and without uplifting it. Once the 
jet has terminated its active phase is rapidly destroyed 
(Reynolds, Heinz & Begelman 2002) and the lobes can 
evolve like buoyant bubbles (see for example Churazov et 
al. 2001 or Quilis, Bower & Balogh 2001). An example 
similar to Perseus A is the Hydra A cluster containing 
the powerful FR I radio source 3C 218. Given the prop- 
erties of the X-Ray gas (a ~ 26 kpc, no ~ 0.06 cm~ 3 , 
T ~ 3.1 keV, David et al. 2001) and a cocoon length 
of ~ 49 kpc we obtain an upper limit on the jet power 
Lj ~ 1.1 x 10 46 ergs -1 . Since there is no indication that 
the gas surrounding the radio lobes is hotter than the am- 
bient cluster gas this limit can be lowered at least by an 
order of magnitude. On the other hand it is possible to 
find some examples of radio sources inside galaxy clusters 
that do not show the presence of X-ray cavities. In our 
scheme these sources could correspond to radio lobes still 
in a strongly overpressured phase. For example the FR II 
radio source 3C 295 is found inside a cluster whose core is 
characterized by a radius a ~ 17.8 kpc, a central density 
no ~ 0.16 cm -3 and a temperature T ~ 3.7 keV (Allen 
et al. 2001). This cluster does not show any cavity in its 
X-ray emission. Given the longitudinal dimension of the 
radio source ~ 17.6 kpc we can estimate a lower limit of 
the jet power Lj ~ 7.1 x 10 45 ergs -1 . 

7. Summary 

In this paper we performed numerical simulations of ax- 
isymmctric supersonic jets propagating in an isothermal 
background atmosphere. In agreement with the results 
presented by Reynolds, Heinz & Begelman (2001), we 
find two distinct and subsequent regimes of interaction 
between the cocoon and the external medium. In the 
first phase of evolution, the overpressured cocoon drives 
a strong shock in the ambient medium, forming a thin, 
hot and compressed shell of shocked material, in the sec- 
ond phase the shock becomes very weak and the shell 
widens, decreasing its density and temperature. The re- 
sulting X-ray morphology in the two phases is different: 
in the strongly overpressured phase, we expect a shell of 
enhanced X-ray emission surrounding the radio emitting 
material, while, in the weak shock phase, we expect a 
deficit of X-ray emission coincident with the radio lobes. 
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We have studied the dependence of the transition between 
these two phases on the physical jet parameters, by a wide 
coverage of the parameter space and by a comparison of 
the results of numerical simulations with analytical mod- 
els. We find that the transition length between the two 
regimes depends essentially only on the jet power scaled 
over a value dependent on the properties of the ambient 
medium. 
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Appendix A: 

As we have discussed, the jet replenishes the cocoon of 
matter and energy, therefore the cocoon expands com- 
pressing the ambient medium and if the expansion speed 
of the cocoon is highly supersonic, it will drive a strong 
shock in the external medium. The first attempts to give 
an analytic description of theses processes was done by 
Begelman & Cioffi (1989) and Ciofii k Blondin (1992). 
To describe the energy input by the jet, they write the 
following energy equation for the extended cocoon 



dE c 
dt 



L 



(AT) 



where E c is the total energy in the extended cocoon and 
Lj is the jet power. Assuming a constant energy input by 
the jet and that all the energy is converted in thermal 
energy, we can then write the average (extended) cocoon 
pressure as 



Pc = 



Vc 



(A.2) 



where V c is the volume of the (cylindrical) cocoon equal to 
7rr 2 Z c , r e is the extended cocoon radius and l c the cocoon 
length assumed to increase as v^t, with Vh given by the 
one-dimensional estimate = ^Jv v-J{\ + \fv), where 
v = jOj/poxt- The lateral expansion of the extended cocoon 
can then be obtained by assuming a strong lateral shock, 
for which we can write 



dr c 
dt 



'7 + 1 



Pcxt 



(A.3) 



where p G xt is the external density. Equation (|A.3|) can be 
integrated giving for the time behavior of r c a dependence 
on t 1 / 2 and for P c a dependence on f -1 . The behavior 
of the cocoon proper can then be obtained by writing its 
energy balance, that will differ from Eq. (|A.1|I by the term 
describing the work done at the contact discontinuity by 
the cocoon on the external medium and will have the form 



dE r 



Li-R 



dV ci 



(A.4) 



dt J dt ' 

where E c \ is now the total energy of the cocoon proper 
and V'ci its volume. This equation can be written in an 
integral form, neglecting the initial value of the volume: 



V C A = 



-J—{l^l f Lf^-DhA . (A.5) 



Inserting in this equation the expression for the pressure 
derived above and assuming that the volume V C i is propor- 
tional to l c rf and the cocoon length behaves in the same 
way as that of the extended cocoon we can derive that the 
cocoon radius r; behaves as t 1 / 2 , in the same way as r c . In 
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this situation wc will then have a thin shell whose relative 
thickness stays constant. This description is well suited for 
strongly overpressured cocoons that drive strong shocks 
in the ambient medium, but the temporal dependence of 
pressure (oc i" 1 ) tells us that a cocoon cannot stay during 
all its evolution in a strongly overpressured regime. In ad- 
dition there might be values of the parameters for which 
the cocoon is not strongly overpressured already at the 
beginning of its evolution. For describing these situations 
we have to extend Begelman & Cioffi (1989) model. We 
need to reconsider the energy balance of the extended co- 
coon, observing that in addition to the energy input by the 
jet, the extended cocoon acquire also the thermal energy 
of the external material that enters the extended cocoon 
during its expansion. We can then write 

Pext 



E, 



L 3 t- 



-dV . 



v c 7 



1 



where the second term on the right hand side represents 
the contribution described above and P ox t is the pressure 
of the external medium. The equation for the pressure, 
instead of Eq. 1A.2(I becomes 

( 7 -l)ij* 



P r = P* 



(A.6) 



PcxtdV . 



(A.7) 



where 

p* = _L 

In the case of uniform external pressure, we have, of 
course, that P* = P cx t- A second extension to the 
Begelman & Cioffi (1989) model is done by considering 
a lateral shock of arbitrary strength and writing the ex- 
pansion speed of the cocoon in a more general way as (see 
e.g. Landau & Lifshitz, 1959) 

/ - — i - i I' \ 

(A. 



dr c 
dt 



7-1 , 7 + 1 Pc 



27 



27 Pext, 

The system of equations (|A.6l - lA~8|) can be solved analyt- 
ically in the case of uniform external pressure assuming 
a cocoon length l c ~ v 1 l 2 v ] t and L } = Tr/2r?pjvf,i.e. the 
kinetic jet power. This approximation is strictly valid for 
jets that are greatly supersonic with respect to their inter- 
nal sound speed, so as to neglect the enthalpy flux term. 
Neglecting its initial value ~ rj (jet radius), the cocoon 
radius then is given by 



c\f + vV " lr^ 1/4 Mc se t = 



Substituting the solution of the pressure in the weakly 
overpressured regime in Eq. I|A.5|I . we obtain that the 
internal radius begins to expand subsonically and tends 
asymptotically to a constant value. 



1 H 



vV 



1 L 



i/V4 r j 
and the cocoon pressure is 



(A.9) 



P. 



7(7 - 1) 



1 



V l/2 M 2 r 2 



i/V 4 M r; 



7(7 - 1) 



I 2 , Vr^T/c 



„3/2 M 4 r 2 „3/4 M 2 

. J J , 



.(A.10) 



